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Abstract 

If the variance of a Gaussian spin-glass Hamiltonian grows like the volume the model 
fulfills the Ghirlanda-Guerra identities in terms of the normalized Hamiltonian co- 
variance. 
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1 Introduction 



In the last decade new ideas and technical methods have been developed in the attempt to 
build the spin glass theory on rigorous mathematical grounds. The most recent example is 
the interpolation strategy [GT] introduced to prove the existence of the thermodynamic 
limit in the Sherrington-Kirkpatrick [SK] model and its use in the proof of the Parisi 
free energy [MPV] of the same model [Gl, Tl] and in the identification of an extended 
variational principle [AiSS]. 

The first basic contribution in the field came with the work by Guerra [G2] on how 
to prove some correlation identities of the Sherrington-Kirkpatrick model that were only 
assumed within the ultrametric structure of the Parisi solution. Those identities were 
later generalized by Ghirlanda- Guerra [GG] and are playing an increasingly important 
role in the mathematical approach to the low temperature spin glass phase (see [B, T2] 
and references therein). The Ghirlanda-Guerra identities are consequence of a very ba- 
sic principle of statistical mechanics i.e. the vanishing of the fluctuation of the energy 
per particle: at increasing volumes the energy per particle approaches a constant with 
respect to the equilibrium measure. Within the non-disordered classical cases all that 
simply implies the finiteness of the specific heat almost everywhere in the temperature 
(see nevertheless the implications in classical mean-field models [CGI]); however in the 
spin glass cases, where the equilibrium quenched state is a properly intertwined composi- 
tion of the Boltzmann-Gibbs and the disorder measures, its consequences are way more 
subtle. The work [G2] led to the identification of the stochastic stability [AC], an invari- 
ance property of the quenched state under a class of suitable perturbations. Subsequently 
stochastic stability was developed and used clarify the relation between the equilibrium 
and the off-equilibrium properties in the spin glass phase [FMPP1, FMPP2]. Recently 
stochastic stability has been classified [C, BCK] and placed on rigorous grounds in [CGi], 
where its relation with the mentioned identities is also discussed. 

In this work we obtain a condition that guarantee the validity of the Ghirlanda-Guerra 
identities: our result states that they hold true whenever the variance of the Hamiltonian 
function grows like the volume. Such a condition is the same that ensure the existence 
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(boundedness) of the thermodynamic limit [CG] and applies to every spin glass model 
studied so far: to the Edwards-Anderson model, to the finite dimensional cases with 
summable or non-summable interactions in the sense of Kanin and Sinai [KS], to the mean 
field cases like SK, p-spin, REM and GREM, up to the general spin glass model of subset 
interaction on which our general condition has been tailored. It is important to stress that 
the identities we prove hold in terms of the normalized Hamiltonian covariance which has 
a different spin expression in each model: for instance in the SK model it coincides with 
the square power of the overlap function, while for the EA model it is the link overlap 
[C2, NS]. The strategy we use to achieve the result relies on very simple methods like the 
bound on martingale sums and classical inequalities. 

The paper is organized in definitions (Sec. 2), results (Sec. 3), proofs (Sec. 4) and is 
concluded with some comments and perspectives (Sec. 5). 

2 Definitions 

We consider a disordered model of Ising configurations o~ n = ±1, n G A C 7L d for some d- 
parallelepiped A of volume |A|. We denote £a the set of all a = {cr n } ne ^, an d |£a| = 2' A '. 
In the sequel the following definitions will be used. 

1. Hamiltonian. 

For every A C 7L d let {-^(oOjo-eEjv be a family of 2' A ' translation invariant (in 
distribution) centered Gaussian random variables defined, in analogy with [RU], 
according to the very general representation 




(2.1) 



where 



(<70 = 0) and the J's are independent Gaussian variables with zero mean 



(2.2) 



Av( J x ) = , 



(2.3) 
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and (translation invariant) variance 

Av(4) = A x 



(2.4) 



2. Covariance matrix. 

C A (a,r) : = Av (H A (a)H A (r)) = ^ A 2 x a x r x . (2.5) 

XcA 

By the Schwartz inequality 

|C A (a,r)| < ^fC^a)^C^7) = ^ A x (2.6) 

XcA 

for all a and r. 

3. Thermodynamic Stability. 

The Hamiltonian (2.1) is thermodynamically stable if it exist a constant c < oo such 
that 

sup ——C A (a, a) = sup — A x < c < oo (2.7) 
AcZ d |A| A C zd IA| - ' 



XcA 



Together with translation invariance a condition like the (2.7) is equivalent to 



X90 

In fact 



XcA xeAX3x 1 1 X30 1 1 

Alternatively, summing over the equivalence classes X of the translation group, the 
(2.7) is equivalent to 

X>| < e. (2.10) 

Thanks to the (2.6) a thermodynamically stable model fulfills the bound 

Ca(<t,t) < c|A| (2.11) 
and has a order 1 normalized covariance 

c A (a,r) := -^C A (a,r) (2.12) 
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4. Random partition function. 



Z(0) := e _/3i/A(CT) . (2.13) 

crtE E A 



5. Random free energy. 

-j3F(j3) := A{(3) := \nZ((3). (2.14) 

6. Random internal energy. 

== y A % A(ff) • (2-15) 

7. Quenched free energy. 

-j3F(J3) := := Av(^)) . (2.16) 

8. R-product random Gibbs-Boltzmann state. 

(3{H A (aW)+-+H A (aW)] 

9. Quenched equilibrium state. 

<-> ■= Av(Q(-)) . (2.18) 

10. Observables. 

For any smooth bounded function G(ca) (without loss of generality we consider 
|Cr| < 1 and no assumption of permutation invariance on G is made) of the covari- 
ance matrix entries we introduce the random (with respect to < — >) R x R matrix 
of elements {qk,i} (called generalized overlap) by the formula 

<G(q)>:= Av(ft(G(c A ))) • (2.19) 

E.g.: G(c A )=c A (a(V (2) )c A (a(V (3) ) 

< 91,292,3 > = Av 



'E CT ( 1 ), CT ( 2 ), CT (3)C A ( ( r( 1 ), ( x( 2 ))c A (a( 2 ) )( x(3)) e ^E;L 1 H A ( 



(2.20) 
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3 Results 



The Ghirlanda-Guerra identities admit several equivalent formulations. They can be ex- 
pressed in terms of factorization properties of the quenched distribution of the generalized 
overlap [GG, B, T2] as well as in terms of expectations of observables. In this work we 
chose the second approach because it allows to distinguish the identities in two classes 
with different physical meaning: the first expresses the regularity with respect to the 
temperature, the second the self-averaging of intensive quantities. 



In relation to the definitions of the previous section it holds the following: 

Theorem 1 The quenched equilibrium state of a thermodynamically stable Hamiltonian 
fulfills, for every observable G and every temperature interval [/5f , /3f ] the following iden- 
tities in the thermodynamic limit 
„p2 R R 

I <J2 G( l^-2RGj2 ( H,R+i + R ( R + 1 ) G( lR + i,R + 2>dp 2 = (3.21) 

J Pi k,l=l 1=1 
k+l 



hi 



R 



< G( lk,R+i > -(-R + l) < Gq R+liR+2 > + < G >< 5i j2 > 



k=l 



dp 2 = (3.22) 



Remark 1 The two previous relations when applied to G(q) = 91,2 combined together lead 
to the well known [MPV, G2]: 

< gi.292,3 > = \ < q 2 > +\ < Q > 2 (3.23) 

1 2 

< 91,293,4 > = 3 < q 2 > +3 < q > 2 (3.24) 

Remark 2 It is straightforward to verify that the condition (2.7) of thermodynamic sta- 
bility holds for all the known spin glass models. Here are a few examples: 

1. The Edwards- Anderson model [EAJ. The nearest neighbor case is defined by = 1 
if X = (n,n') and \n — n'\ = 1. The condition (2.7) is verified by c = d. 
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2. More generally one consider still a two body interaction with A x = \n — n'\ 2da . The 
regime a > 1 comes from a summable interaction. The condition (2.7) is verified 
byc = (2a - l)~ d for all a > 1/2 thus including also the non summable case [KS]. 

3. The SK model [SK]. Although it is not a finite dimensional model it may still be 
embedded in 7L, with A\ = unless \X\ = 2 and A it j = N^ 1 with N = |A|. It 
obviously fulfills condition (2.7) with c = 1 . 

4- The p-spin. Analogously as above Ax = if \X\ ^ p and A x = 1/N P otherwise. It 
is thermodynamically stable with c = 1 

5. The REM [D] and GREM [EG] models. Although they have not been defined as 
spin models their discrete nature allows to associate to them a spin Hamiltonian. 
For instance it is easy to prove that the REM is represented by by the Hamilto- 
nian (2.1) with A 2 X = N2~ N which satisfies the condition (2.7) with c = 1, see 
also [Bj. The same argument holds for the GREM [CDGG, CG2] which is again 
thermodynamically stable with c = 1. 

Remark 3 The relevance of the identities is evident considering that they reduce the 
degrees of freedom a priori carried by each spin glass model. In the mean field case 
for instance the method led to the rigorous proof of a property called replica equivalence 
[Pa, C2] which can viewed as an ansatz generalizing the ultrametric one. The purely 
ultrametric identities (still lacking a rigorous mathematical derivation) which are built in 
the Parisi solution of the SK model are not contained in the Ghirlanda-Guerra ones. 

Remark 4 It would be interesting to establish, or disprove, the same identities in a 
stronger sense, i.e. everywhere in the temperature. One of the limit of the method used to 
achieve our results is that it is intrinsically restricted to hold in [3-average, i.e. in every 
interval excluding at most isolated singularities. It is still an open question if, in the spin 
glass phase, there are similar singularities or if the identities hold just everywhere. The 
only existing results are evidences of numerical nature of the validity of those identities 
everywhere [MPRRZ, CGi2j. 
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4 Proof 



The statements (3.21) and (3.22) are proved respectively in the lemmata of subsections 4.1 
and 4.2. The proof uses only elementary methods like martingale differences and classical 
inequalities. Let h(a) = | A^i^o") be the Hamiltonian per particle. We consider the 
quantity 

R 

{< % (0 ) G > - < h(a^) >< G >} = AiG + A 2 G (4.25) 

i=i 

where 

R 

A 1 G = { Av {tt[h(a^) G] - n[h(<r®)]n[G\)) } (4.26) 



i=i 

R 



A 2 G = J2 ( Av (n[h(<r®)]n[G\) - Av (n[h(<r®)]) Av (n[G\)} (4.27) 



1=1 



4.1 Stochastic Stability Bounds, vanishing of AiG 

We follow the method of stochastic stability as developed in [CGi]. 

Lemma 4.1 For every bounded observable G, see definition (10), we have that for every 
interval [/^i,/^] in the thermodynamic limit 

/ A ± Gdp = (4.28) 

Proof. We observe that deriving < G > with respect to the temperature 

R 

= \M E ( Av ( n [M^ (0 ) G] - m° {l) )MG\)) } (4-29) 
Integrating in d(3 we obtain thanks to (4.26) 

r\ iGdp= ^G>( M -<G>w (4 30) 

hi l A l 

Remembering the assumption on boundedness of function G this proves the lemma. □ 

Remark 5 The previous lemma is related to a general property of disordered systems 
which is known as stochastic stability (see [AC, CGi]). It says that the equilibrium state 
in a spin glass model is invariant under a suitable class of perturbation in all temperature 
intervals of continuity. 



Lemma 4.2 The following expression holds: 



AiG = -p < G 



R 



R 



J2 qi,k-2Rj2 ( li,R+i + R ( R + l ) 

<lR+l,R+2 



k,l=l 
L k^l 



1=1 



> 



(4.31) 



Proof. 

For each replica / (1 < I < R), we evaluate separately the two terms in the right side 
of Eq. (4.26) by using the integration by parts (generalized Wick formula) for correlated 
Gaussian random variables, x±, X2, ■ ■ ■ , x n 

Av (a* iK*!, ...,*„)) = |>v Av f W(xu~.,x n ) \ (4 32) 

j = l \ 3 J 

It is convenient to denote byp (R) the Gibbs-Boltzmann weight of R copies of the deformed 
system 



p(R) 



R 



so that we have 



(4.33) 



(4.34) 



We obtain 



Av (n(h(a®)G)) 



y^Avj GH A (*V)p(R)\ 



av y: e gc ^ {i) ^ 



aW,...,a( r) T 
R 



dp(R) 
dH A (r) 



Y < Gqi, k > -R < Gq liR+1 > 



k=l 



(4.35) 
(4.36) 
(4.37) 



where in (4.36) we made use of the integration by parts formula and (4.37) is obtained 
by (4.34). Analogously, the other term reads 



Av(n(M<rW))n(G)) = ^rAvfx; Yl GH A (*W)p(R+l)\ 

1 1 \ aW r(D,...,r(«) / 



(4.38) 
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-f3 J2 <G( 1^r + i>-(R+1)<G QR+l,R+2 > 




(4.39) 



_fc=i 



(4.40) 

Inserting the (4.37) and (4.40) in Eq. (4.26) we finally obtain the expression (4.31). □ 
4.2 Martingale Bounds, vanishing of A2G 

The method of the martingale differences to prove the self averaging of the free energy, or 
in general to bound the fluctuations of extensive quantity has been applied in the context 
of spin glasses in [PS] for the SK case and in [WA] in the case of finite dimensional models. 
Our formulation applies to both cases and extends the previous results. For instance our 
method includes the non summable interactions in finite dimensions [KS] and the p-spin 
mean field model as well as the REM [D] and GREM [DG] models. 

Lemma 4.3 The free energy is a self averaging quantity, i.e. it exist a positive function 
c(/3) such that 



Proof. For an assigned volume A we enumerate by the index k the interacting subsets 
X from 1 to and considering the random partition function (2.13) we define 



where the symbol Av< fc denotes the Gaussian integration performed only on the first k 
random variables Jx- Clearly A = A(f3) and A Nh = A(f3). Introducing the quantity 



V(A) = Av (A 2 ) - Av(A) 2 < c(/3)|A| 



(4.41) 



(4.42) 



= A k -A, 



k+l j 



(4.43) 



it holds 



JVa-1 




(4.44) 
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and 

V(A) = E Av (^) + 2^Av(W fc ,) . 

k k>k' 

First we observe that the second sum is zero, being zero each of its terms. In fact 
Av(* fc *fcO = Av (Av< fc (* fc *fcO) = Av(* fc Av< fc (*fcO) 

and 

Av< fc (* fc /) = Av< k (A k , - A k , +1 ) = 

thanks to the property 

Av< k (A k ,) = A k V k>k'. 
We introduce now the interpolated Hamiltonian 

N A 

h a(°) = -J2 Jitiai 
1=1 

with 

t, if I = k + 1, 



*i = 

and define the quantity 



1, otherwise, 



By the fundamental theorem of calculus 

A k = A k {0) + B k 

with 



B k 

We observe 



f 1 dA k (t) , /* ^ A 



Av(tf£) = Av(L4 fc - A fc +i] 2 ) = Av(Av fc+ i[A fc - Av k+1 A k } 2 ) 

= Av{Av k+1 (A 2 k )-[Av k+1 (A k )] 2 ) 



11 



Since A k and B k differ by a constant with respect to Avfc+i (integration with respect to 
the (k + l)-th Gaussian) we have that their variance is the same: 

Av(*y = Av(Av fc+1 (^)-[Av fc+1 (A,)] 2 ) = Av(Av fc+1 (^)-[Av fc+1 (fi fc )] 2 ) . 

(4.55) 

We will estimate separately the two terms Av (Av fe+ i(5|)) and Av ([Av k+ i(B k )] 2 ). By a 
simple integration by parts (4.32) on J k+ i we obtain 

Av k+1 (B k ) = P 2 *l+if Av< k+1 [l -<j*(a k+1 )]tdt < \p 2 A 2 k+l (4.56) 

which implies 

< Av ([Av k+1 (B k )] 2 ) < \p A Ai +1 . (4.57) 

Analogously we have 

Av fc+ i(B^) = Av fc+1 / / A\< k (uj t (J k+1 a k+1 ))Av< k (uj s (J k+1 a k+1 ))stdsdt (4.58) 

Jo Jo 

Applying twice the integration by parts (4.32) we get 

Av k+1 (Bl) = (3 2 A 2 k+1 Av k+1 [ [ Av< k (cu t ((T k+1 ))Av< k (uj s (a k+1 ))stdsdt + (4.59) 

Jo Jo 

- 2p 4 A 4 k+1 Av k+1 / Av< k (uj t (a k+1 )[l - uj 2 (a k+1 )])Av< k (uj s (a k+1 ))st 3 dsdt + 

Jo Jo 

- 2^A\ +l Av k+l [ [ Av< k (u s ((r k+1 )[l-uj 2 s (a k+1 )])Av< k (u; t (a k+1 ))shdsdt + 

Jo Jo 

+ 2p 4 A 4 k+1 Av k+1 [ [ Av< k [l-u 2 (a k+1 )]Av< k [l-ij 2 (a k+1 )}s 2 t 2 dsdt 
Jo Jo 



Putting together the (4.57) and the (4.59) we find 

Av (*y < \p 2 A 2 k+1 + f^ A At +1 (4.60) 

V(A) = E Av (*» < E + |^ • (4.61) 

k XcA 

By the assumption of thermodynamic stability with the formulation (2.10) and using the 
inequality J2x A 4 ~ < (J2x ^x) 2 we obtain 

1 ^ 

V(A) < |A|( i /3 2 c+-/3 4 c 2 ) (4.62) 
12 



which fulfills (4.41) with c(P) = \f3 2 c + f/9 4 c 2 . □ 

Lemma 4.4 The internal energy is self averaging almost everywhere f3, i.e. defining 
u = U/\A\ and V(u) = Av (u 2 ) — Av (u) 2 it holds in the thermodynamic limit 

V(u)d(3 -> (4.63) 



/ 



•Pi 
Proof. 

The result is obtained in two steps which use general theorems of measure theory. First 
from lemma 4.3 we obtain the convergence to zero almost everywhere (in /3) of the variance 
of the internal energy, then thanks to a bound on the variance of the internal energy we 
apply the Lebesgue dominated convergence theorem which gives the lemma statement. 
The sequence of convex functions A((3)/\A\ converges a.e. (in J) to the limiting value 
a(/3) of its average [GT, CG] and the convergence is self averaging in the sense of lemma 
4.3. By general convexity arguments [RU] it follows that the sequence of the derivatives 
A'((3)/\A\ converges to u(f3) = a'(/3) almost everywhere in (3 and also that the convergence 
is self averaging. In fact the vanishing of the variance of a sequence of convex functions 
is inherited, in all points in which the derivative exists (which is almost everywhere for a 
convex function), to the sequence of its derivatives (see [S, OTW]). From lemma 4.3 we 
have then 

V(u) -> (3 - a.e. (4.64) 

In order to obtain the convergence in /5-average we use the Lebesgue dominated con- 
vergence theorem. In fact we prove that the sequence of variances of u is uniformly 
bounded (in every interval [Pi, P2]) by an integrable function of (5. A lengthy but simple 
computation which uses again integration by parts gives 



Av(W) = Av(j2 J x^x)) = ^T/^l-Av^Vx))] <!3\A\c (4. 

VxcA / XcA 

Av{U 2 ) = Av I ^ J x J y lu(o- x )uj(o-y) ) = 

\X,YCA J 



65) 



(4.66) 

f3 2 A 2 x A 2 Y Av [1 - uo 2 {a x ) - uo 2 {a Y ) + Quo 2 {a x )uo 2 {a Y ) + 

X,YCA 

Qlu(o- x )lu(o-y)lu(o- X (Ty) +LU 2 (a x a Y )] < U(3 2 \A\ 2 c 2 (4.67) 
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from which 



From this follows (4.63). 



V(u) < 15p 2 c 2 



(4.68) 
□ 



Lemma 4.5 For every bounded observable G, see definition (10), we have that for every 
interval P2] in the thermodynamic limit 

"02 



[ A 2 Gd(3 = 



(4.69) 



Proof. 

Thanks to the Schwartz inequality 

A 2 G = Av (uG - Av (u) Av (G)) — Av ([u - Av (u)] [G - Av (G)]) — (4.70) 



< JAv([u- Av (u)] 2 ) J Av ([G - Av(G)] 2 ) < y/2y/V(u) 



(A 2 G) 2 < 2V{u) 



02 



A 2 Gdp 



< 



I f02 / f02 

J J (AsGO^v/^A < VW2-Pi)J J V(u)df3 



(4.71) 


(4.72) 
□ 



Lemma 4.6 The following expression holds: 



A 2 G = -PR 



R 



< Gq l>R+1 > ~(R + l) < Gq R+hR+2 > + < G><q h2 > 



k=i 



(4.73) 



Proof. In order to obtain the A 2 G we are left with the explicit evaluation of the other 
term in (4.27) which simply gives 

Av(Q(% (0 ))) Av(Q(G)) = 4jAv( £ H A (a^)p A (l)) <G> 

= av(x;e c a(^7) 



dpA (1) 
^a(7) 

= -P < G > [< g M > - < g 1>2 >] 



< G > 



(4.74) 
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Inserting the (4.40) and (4.74) in Eq. (4.27) we obtain the (4.73). 



□ 
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